The influence of chemical composition on models of Type la supernovae 
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Type la supernovae are bright stellar explosions distinguished by standardizable light curves that 
allow for their use as distance indicators for cosmological studies. Despite the highly successful use 
of these events in this capacity, many fundamental questions remain. Contemporary research inves- 
tigates how properties of the progenitor system that follow from the host galaxy such as composition 
and age influence the brightness of an event with the goal of better understanding and assessing 
the intrinsic scatter in the brightness. We provide an overview of these supernovae and proposed 
progenitor systems, all of which involve one or more compact stars known as white dwarfs. We 
describe contemporary research investigating how the composition and structure of the progenitor 
white dwarf systematically influences the explosion outcome assuming the progenitor is a single 
white dwarf that has gained mass from a companion. We present results illustrating some of these 
systematic effects from our research. 

PACS numbers: 97.60. Bw,26.50.+x,97.20.Rp,26.30.-k,95.30.Lz 



I. INTRODUCTION 

A. Historical Overview of Supernovae 

"New" or "guest" stars have occasionally appeared in 
the sky since time immemorial. Ancient observations in- 
clude solar system objects such as comets and stellar ex- 
plosions that appear suddenly as a newly visible star. 
The earliest record of an event associated with a stel- 
lar explosion is found the The Book of the Later Han 
which describes the appearance of a "guest star" in the 
year 185 [ILQ. Modern study of stellar explosions began 
with Tycho Brahe's naked-eye observations of the event 
we now call Supernova 1572, and the name "nova" ap- 
plied to the phenomena of stellar explosions follows from 
his book De Stella Nova ("On the New Star"). Tycho's 
observations were of such fidelity that later astronomers 
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were able to reconstruct the light curve and color evolu- 
tion of this event for comparison and classification 0, 0] ■ 

By the early twentieth century, technological advances 
in astronomy led to recognition of the vast distances to 
other galaxies, implying that novae observed in other 
galaxies were in fact extremely bright. This advance, 
together with advances in theory, resulted in the classi- 
fication of "super-novae" as bright explosions associated 
with the violent death of a star (Ref. H; see also Refs. @ 
and0 for discussion of pioneering work). Observations of 
a supernova typically produce spectra and a light curve, 
with spectra taken during the course of an event pro- 
viding information identifying the constituent elements 
and the light curve plotting the intensity of light vs. 
time. As originally proposed by Minkowski Q, super- 
novae are classified by properties of their spectra and 
light curves, with the principal distinction arising be- 
tween events without evidence of hydrogen in their spec- 
tra (Type I) and events with evidence of hydrogen in 
their spectra (Type II). Type I supernovae were further 
divided into Types la, lb, and eventually Ic 0-[T3- The 
Type la sub-classification depends on the observation of 
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a specific Si line not seen in Types lb and Ic. Type 

lb and Ic are distinguished by Type lb spectra exhibiting 
He features while Type Ic spectra do not. 

Given the extreme brightness of supcrnovae, only two 
possible energy sources are thought to account for these 
phenomena: the release of gravitational potential energy 
or the release of nuclear binding energy. Events in all 
of the observational classifications, with the exception of 
Type la, are understood to follow from the gravitational 
collapses of a massive star that has exhausted its nu- 
clear fuel [l3| . The remnant in this type of explosion is 
a neutron star or black hole, of which many have been 
observed. Compact remnants have never been observed 
from Type la supernovae. While the classification fol- 
lows from the spectral signature, Type la supernovae are 
now understood to be the result of a thermonuclear ex- 
plosion consuming roughly one and a half solar masses 
of degenerate stellar material composed principally of C 
and O @. The progenitor systems of these events, how- 
ever, remain the subject of debate. Most proposed Type 
la progenitor systems follow from the suggestion of Hoyle 
and Fowler [14| of a thermonuclear runaway occurring in 
the core of a star supported by electron degeneracy. Var- 
ious settings in which suitable conditions for the explo- 
sive burning of degenerate stellar material have been ex- 
plored, and scenarios typically involve one or more white 
dwarf stars, as we discuss below. 



high (z ~ 1) rcdshifts, the systematics of this evolution 
must be accounted for to allow precision cosmological 
measurements. Avoiding relying on empirical relation- 
ships is one important goal of the theoretical study of 
Type la supernova explosions. 

While the yield of 56 Ni explains the first-order vari- 
ations in the light curve, current research also aims to 
understand higher-order effects and the physics behind 
the Phillips relation. Observations report that Type la 
supernovae appear to have an intrinsic scatter of a few 
tenths of a magnitude after calibration, forcing a min- 
imum uncertainty in any distances measured by using 
Type la supernovae as standardizable candles [3, [25j . 
An important goal of theoretical research into Type la 
supernovae, from the standpoint of cosmology, is to un- 
derstand the sources of scatter and to identify potential 
systematic biases by studying the effects of various prop- 
erties on the mechanism and nucleosynthetic yield of the 
supernova. The surrounding stellar population and vari- 
ous properties of the progenitor such as its composition, 
zero-age main sequence mass, thermodynamic state prior 
to ignition, and cooling and accretion history are known 
to affect the lightcurves of Type la supcrnovae; the role 
of these "secondary" parameters is the subject of consid- 
erable study (e.g., Refs. [26|-[28j). Additionally, many of 
these effects may be interconnected in complex ways |29l — 




B. Observations of Type la Supernovae 



C. Trends in Observations of Type la Supernovae 



The release of nuclear binding energy via explosive 
burning of degenerate C and O provides more than 
enough kinetic energy to unbind a white dwarf star, and 
the expansion velocities of freshly synthesized elements 
typically reach on the order of 10, 000 kms^ 1 . The peak 
brightness of a Type la supernova is set not by the ex- 
plosion energy, but by the synthesis in the explosion of 
radioactive 56 Ni, which decays via the chain 56 Ni to 56 Co 
to 56 Fe releasing photons that power the observed light 
curve 

The vast majority of Type la supcrnovae obey a cor- 
relation in which the peak brightness is positively corre- 
lated with the timescale over which the lightcurve decays 
from its maximum [l9l | . This "brighter is broader" trend 
is known as the Phillips relation [2(| and it allows the 
peak brightnesses to be calibrated so that these events 
may be treated as "standard candles" for determining 
distances. The resulting parameter describing the scaling 
is known as Amis^B), which is the change in the appar- 
ent B-band magnitude from peak brightness to 15 days 
later. The correlation is understood physically as stem- 
ming from having both the luminosity and opacity b eing 
set by the mass of 56 Ni synthesized in the explosion [21- 
[23j . This relation has been exploited to make Type la 
supcrnovae the premier distance indicators for cosmolog- 
ical studies. We note that because the galaxy /stellar 
population is known to have evolved significantly since 



The interest in Type la supernovae that follows from 
the successful use of these events as distance indica- 
tors in studies revealing the acceleration of the Uni- 
verse [H, [HI has led to modern observational campaigns 
designed to ex plore sy stematic effects on the brightness 
of these events [35l - l37| . Observations explore correlations 
between observed properties of an event such as peak 
brightness and properties of the host galaxy such as its 
composition, age, and mass. Many of these galactic prop- 
erties are correlated, which makes the task of determining 
the underlying physical reason for correlations between 
properties of supernovae and their host galaxies difficult. 
Compounding this difficulty is the fact that detailed ob- 
servations of galaxies are not available at high redshifts, 
and parameters (typically galaxy mass and star forma- 
tion rate) are therefore inferred from galactic models that 
reproduce the observed spectral energy distribution. 

Hamuy et al. j^Ij reported that the peak brightness 
correlates with the morphological type of the host galaxy, 
with the most luminous supcrnovae tending to be hosted 
by late-type, spiral galaxies. While the trend is striking, 
supernovae within the same morphological type do not 
necessarily share the same physical environment. Several 
physical properties of galaxies tend to correlate with the 
morphological type. Early-type, elliptical galaxies are 
thought to form through galaxy mergers and are more 
massive, contain older stellar populations, and have lit- 
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tie star formation. The composition of a galaxy depends 
largely on the proportion of material that has previ- 
ously been processed in stars (i.e. elements other than H 
and He, which are collectively referred to as "metals"). 
"Metallicity," the relative abundance of these elements, 
in the gas phase is correlated with the galaxy mass as 
galaxies with deeper gravitational potential wells tend 
to more effectively retain metals [39|. Therefore, early- 
type galaxies also tend to have higher metallicities. Late- 
type, spiral galaxies are less massive, are actively form- 
ing stars, and hence contain younger stellar populations. 
Likely some combination of these properties bias the pro- 
genitors of Type la supernovae to produce the observed 
correlation with galaxy morphology. 

Despite challenges, Gallagher et al. [Z(| observed a 
slight dependence of peak brightness on metallicity with 
dimmer events in metal-rich galaxies. A conclusive trend, 
however, has proven elusive. Due to the weak correla- 
tion with metallicity, progenitor age is suspected to be 
primarily responsible for variations in the peak bright- 
ness [40l - [43l ]: however, challenges remain in drawing this 
conclusion. First, the present gas-phase metallicity of the 
host galaxy may not be representative of the metallicity 
of the gas the progenitor star formed from, particularly 
for older stellar populations. Second, the determination 
of the age of a stellar population from the integrated 
light of a galaxy is "de gen erate" with the determination 
of the metallicity jiil l45j|. That is, both effects tend 
to redden the light from a galaxy, so it is difficult to 
distinguish the true cause of the reddening. For obser- 
vations of high-rcdshift galaxies, the metallicity must be 
inferred from the galaxy mass, precluding the possibility 
of directly determining its effect on the supernova light 
curve [43]. When correlations with the inferred stellar 
population age are considered, dimmer supernovae are 
found in older populations. 

Metallicity could also be responsible for variations 
in the peak brightness of calibrated light curves. At 
present, observers do not use the Phillips relation [2(| but 
employ much more sophisticated methods to calibrate 
light curves and thereby measure extragalactic distances. 
These methods construct calibrated light curves in mul- 
tiple pass-bands based on templates from nearby events 
and include corrections for extinction, time dilation, and 
the K correction for shifting of spectra due to cosmo- 
logical redshift. Examples include the Spectral Adap- 
tive Light Curve Template SALT [H, 53, th e Multicolor 
Light-Curve Shapes method MLCS2k2, [48|,|49|, and the 
(Python-based) supernova light curve fitter SNooPy [5(| ■ 

Any metallicity dependence reported for calibrated 
light curves is representative of systematic effects that 
have not been taken into account in the calibration pro- 
cedure. Gallagher et al. [4(1 El| report a dependence 
in the peak brightness in the calibrated light curves with 
metallicity, suggesting the possibility of systematic effects 
with redshift that would have implications on cosmologi- 
cal parameters and the equation of state of Dark Energy 
Howell et al. [43|, however, using a different calibration 



method, find no such dependence on metallicity. 

Theoretically, the presence of metals in a progenitor 
white dwarf influences the outcome of the explosion by 
changing the path of nuclear burning, which influences 
the amount of 56 Ni synthesized in an event [5l(. Be- 
cause the decay of 56 Ni powers the light curve, metallic- 
ity can therefore directly influence the brightness of an 
event. Metallicity can also influence the outcome of an 
explosion in other ways, including changes in the struc- 
ture of the white dwarf that result from the influence of 
metallicity on stellar evolution, sedimentation within the 
white dwarf, the nuclear flame speed, additional sources 
of opacity, and properties of the thermonuclear burning 
such as the ignition density (described below) [2^, [52j . 

In addition to properties of the light curves, Type la 
supernovae rates may also yield clues to the progenitors 
of these explosions. Unlike other types of supernovae, 
Type la have a non-zero rate in older populations, which 
also points to a white dwarf progenitor. Some results in- 
dicate that the dependence of the la rate on stellar age 
is best fit by a bimodal distribution having a prompt 
component less than 1 Gyr after star formation and a 
tardy component several Gyr later [53], l54j . with the 
prompt component appearing brighter on average than 
the tardy component. A bimodal distribution with delay 
time or mean stellar age would imply two distinct pro- 
genitor formation channels. Other studies only indicate a 
correlation between the delay time and the brightness of 
Type la supernovae favoring a single formation channel, 
with dimmer events occurring at longer delay times [4lT - 

Perhaps the most challenging aspect of observing su- 
pernovae is accounting for dust. The process of calibrat- 
ing light curves and inferring luminosity distances relies 
on using the width-luminosity relation and color correc- 
tions, and the manner in which color is linked to dust is 
considered controversial (56|. Color variations observed 
between supernovae correlate with the dust content of the 
host galaxy, but it is also possible that some of the color 
variation is intrinsic to properties of the explosion. Ob- 
servationally disentangling intrinsic variability from ef- 
fects of dust is the subject of active research [56( . 

Observations find a surprising range in the variation 
of the intrinsic brightness of these events, including very 
bright events that suggest more material burned than 
could originate from a single white dwarf [57l - l60j . The 
logic of observations suggesting that more material burns 
than could come from a single white dwarf is based on the 
fact that a white dwarf is supported by electron degener- 
acy and is therefore subject to a maximum mass, known 
as the Chandrasekhar limit, beyond which it will col- 
lapse to a neutron star [6lT - [63l | . A bright observation that 
suggests more 56 Ni synthesized than the Chandrasekhar 
maximum mass therefore implies more than one white 
dwarf was involved. Accordingly, contemporary research 
explores the efficacy of several progenitor systems with 
one or more white dwarfs. 
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D. Proposed Progenitors Systems of Type la 
Supernovae 

We briefly introduce some proposed progenitor systems 
of Type la supernovae and outline contemporary research 
in each. The accepted result for a Type la supernova is 
an event that produces about O.6Af of 56 Ni, the decay 
of which produces the light curve, and most models posit 
the thermonuclear burning of a mix of C and O under de- 
generate conditions (see Refs. [§4|, [65J for examples of al- 
ternative models.) Reviews of progenitor models may be 
found in Hillebrandt and Niemeyer [f|, Livio [6a |. Wang 
and Han While we note that events produce O.6M 
of 56 Ni, there is of course scatter in the 56 Ni yield (and 
brightness) of events and sub-classifications, e.g. Branch- 
normal, 1998bg-, 1991T-, and 2002cx-like events (gUzil- 

The association of Type la supernovae with white 
dwarfs began with the work of Hoyle and Fowler [l4[ who 
found that the thermonuclear incineration of a mixture of 
C and O under degenerate conditions could explain Type 
I supernovae. Hoyle and Fowler theorized a mass 
approaching the Chandrasekhar limit would be required, 
but the question remained as to where an appropriate 
massive stellar core might be found. A suitable progeni- 
tor system must satisfy two requirements. First, the sys- 
tem must have enough mass to produce the amount of 
intermediate-mass and Fe-group elements, e.g. 56 Ni, in- 
ferred from observations, a small fraction of which imply 
more synthesized mass than could come from a single 
white dwarf. Second, while massive white dwarfs may 
be composed of mixtures of C and O or O, Ne, and Mg, 
the progenitors of Type la supernovae must be composed 
principally of C and O because the temperatures required 
to ignite O or Ne cannot be achieved prior to collapse (see 
Ref. [72j and references therein). Accordingly, proposed 
progenitor systems typically involve either multiple C-0 
white dwarfs or mass transferred from a companion onto 
a C-0 white dwarf. 

Many questions concerning progenitor channels remain 
and are actively being explored as described below. Con- 
straints on progenitor channels are also under study from 
the perspective of population synthesis (see for exam- 
ple Refs. (73l - [7r| ). We note that many questions remain 
concerning the fundamentals of thermonuclear burning 
in white dwarf material [77j . and, as mentioned above, 
disentangling the primary from secondary parameters in 
the observed light curve is a challenge. These many issues 
influence the question of Type la progenitor channels. 



1. A Single Massive White Dwarf 

This progenitor channel is referred to as a "single de- 
generate" due to the basic requirement that the white 
dwarf accretes material from a non-degenerate compan- 
ion such as a red giant, helium, or main sequence star 
[7§ ] (see also Refs. |79l - [82l |). In fact, the single degenerate 
channel may be comprised of more than one evolution- 



ary pathway, and several are actively being investigated. 
While no particular channel has direct observational ev- 
idence supporting it, some evolutionary pathways are 
scrutinized more heavily due to the lack of expected ob- 
servational features in the light curve and spectra. In 
particular, a red giant companion to an exploding white 
dwarf is expected to produce a detectable feature in the 
UV when the blast wave interacts with the red giant gas 
[83^ To date, observations of this type of interaction are 
lacking, which rules out the possibility that the red giant 
channel is the only evolutionary path. It may still be 
possible, though, that this pathway contributes a small 
percentage of the total rate of Type la supernovae. 

Some specific observations of Branch-normal events 
provide upper limits of the mass of the companion of 
those particular systems. For example, the recent SN 
201 lfe is the closest Type la supernova to explode in 
modern times allowing the earliest time observations as 
well as detailed pre-cxplosion images that provide an up- 
per limit to the brightness of the companion star [HI, [H| . 
These observations have ruled out luminous red giants 
and almost all helium stars as the companion to this par- 
ticular system [86| . Wheeler [87j has suggested M-dwarfs 
as a viable donor star due to the expected rates at which 
WD /M-dwarf binaries are produced; however, further re- 
search is needed to explore the unique accretion process 
that is assisted by the magnetic properties of both the 
M-dwarf and white dwarf. The evidence in PTF llkx 
[111 of nova shells predating the supernova explosion im- 
plies that at least some fraction of fairly normal Type la 
supernovae must originate from single degenerate type 
systems. 

The single degenerate W7 model by Nomoto et al. [52j 
that used a parameterized burning velocity offered good 
agreement with observations and strongly motivated fur- 
ther study of the single degenerate progenitor. The key 
to the success of W7 follows from the fact that the den- 
sity at which degenerate stellar material burns largely 
determines the composition of the ash. At high densi- 
ties (p > 10 7 gcm~ 3 ), nuclear fuel reacts by fusion all 
the way to the Fe group on the timescale of the burning 
front propagating though the white dwarf. At lower den- 
sities, the burning is incomplete and stops at the Si group 
or even after C burning. The parameterized reaction of 
W7 allowed the star to expand prior to the completion of 
burning, which produced a stratified remnant that agreed 
well with observations (see Fig. 4 of Ref. [89(). 

As the success of the W7 model demonstrates, the na- 
ture of the burning in the single degenerate progenitor 
scenario is critical to the explosion outcome. Arnett (90j 
investigated C detonations, a case in which burning pro- 
ceeds as a supersonic front, in the C-0 cores of massive 
stars. The result was that the entire core was incinerated 
to the Fe group, which is now understood as producing 
too much Fe. The alternative to a detonation is a de- 
flagration, a subsonic burning front that, in the case of 
white dwarf material, proceeds by the conductive prop- 
agation of heat. Factors such as the interaction of the 
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flame with turbulence influence the burning rate, but af- 
ter considerable study, the consensus is that a pure de- 
flagration will not produce a normal Type la supernova 
because the relatively slow burning of the deflagration ex- 
pands the star too much to match observed abundances 
(see Ref. [§lj] and references therein). 

The favored single degenerate explosion model is that 
of a period of deflagration (subsonic burning) that al- 
lows the star to expand followed by a detonation (su- 
personic burning) that rapidly incinerates the expanded 
star. This delayed detonation paradigm was introduced 
by Khokhlov J92||, and simulations of this class of models, 
like Nomoto's W7, demonstrate good agreement with ob- 
servations IHjHH and have been widely accepted. Many 
variations on this theme exist, however, including pul- 
sational detonations in which the white dwarf expands 
and recollapses, subsequently igniting a detonation (9^ - 
gravitationally confined detonation (GCD) in which 
compressional heating ignites a detonation in material 
flowing across the surface of the white dwarf [53, HH (see 
also Refs. [99l . llOOf ). and models in which a deflagration 
transitions to a deto nation w hen the right local condi- 
tions arc met 

Questions remain about the details of flame ignition: 
does a localized region initiate explosive burning or is it 
distributed? If localized regions runaway, how many are 
there and how are they distributed? Due to the highly 
nonlinear evolution of the deflagration phase, the details 
of the ignition process are critically important for deter- 
mining the outcome of an explosion. At one extreme, 
single, off-center ignitions likely lead to an asymmetric 
deflagration phase, like that described in the GCD sce- 
nario, while distributed central ignitions lead to fairly 
symm etric deflagrations. A recent study by Maeda et al. 
jl07| | suggested that distributed, off-center ignitions pro- 
duce viewing angle effects that influence the peak bright- 
ness, and may be largely responsible for the observed 
variation. 



2. Merging White Dwarfs 

The idea of merging white dwarfs as the progenitors of 
Type la supernovae has been around for some time with 
investigations into near- an d super-Chandrasekhar mass 
coalesced objects 108l4lll |. The simplest picture of the 
merging white dwarf scenario has tidal effects disrupting 
the secondary, which forms a disk from which the primary 
subsequently accretes C and O. One of the most com- 
monly cited perils of the merging white dwarf model is 
that the C will ignite at the edge of the coalesced object, 
and a flame will burn inward, conv erting the C-0 white 
dwarf into an O-Ne-Mg white dwarf [112l ll 131 ] . Instead of 
becoming a Type la supernova, further accretion leads to 
the collapse of the white dwarf into a neut ron s tar, a sce- 
nario termed "accretion-induced collapse" 



lapse can be avoided if t he accreti on rate is low enough 
(< few xlCT 6 M yr" 1 [ill Ell) that the heating at 
the surface of the primary does not ignite the C. Models 
of violent mergers of white dwarfs (described below) also 
can avoid an accretion-induced collapse by initiat ing a 
detonation of the sub-Chandrasekhar mass primary [116[. 
Accordingly, there is great interest in demonstrating that 
the detailed accretion and burning processes provide a 
mechanism to avoid accretion-induced collapse. 

There are a number of studies reported in the literature 
that have advanced the state of understanding of white 
dwarf mergers that we highlight here. Other st udies in- 
vestigated collisions of white dwarfs |H7l - ll20l |. These 
events may produce Type la supernovae, but because the 
likely locations for such collisions are in the dense cores 
of globular clusters, these events will be rare, perh aps a s 
few as 10-100 per year in the local (z < 1) universe [l!8l |. 
There is also a core-degenerate merger scenario in which 
a white dwarf merges wit h the core of an asymptotic gi- 
ant branch star (see Ref. jl2l| and references therein). 

Benz et al. [l22| performed some of the earliest simula- 
tions of merging white dwarfs and found that the merger 
remnant is best described as a hot, rapidly rotating thick 
disk consisting of the tidally-disrupted secondary sur- 
rounding the original primary white dwarf. These sim- 
ulations confirmed the basic idea that the inspiral of 
the white dwarfs driven by gravitational radiation can 
cause the ultimate merger of the system, and empha- 
sized that angular momentum transport is the primary 
physical mechanism that shapes the outcome. 

Yoon et al. [123{ modeled the merger of 0.9 Af Q and 
0.6 Mq white dwarfs and argued that accurate modeling 
of the merger event is critical to determining the details 
of the coalesced remnant. Their simulation indicated the 
less massive star was disrupted, and the result was de- 
scribed as "a differentially rotating single C-0 star con- 
sisting of a slowly rotating cold core and a rapidly ro- 
tating hot extended envelope surrounded by a Kcplerian 
disk" instead of the white dwarf and thick disk found 
in earlier studies. Based on mapping the result of the 
three-dimensional simulation to one dimension to study 
burning, they argued that this configuration does not 
give rise to ignition at the edge of the star and hence the 
formation of an O-Ne-Mg white dwarf. Instead, central 
ignition of the C-0 remnant can occur. 

An extensive study was made by Loren-Aguilar et al. 
124j . who considered different systems with varying 



masses. As with the Yoon et al. case, they found that the 
systems avoid igniting the C explosively at the instant of 
merger. They found a coalesced remnant configuration 
similar to that of Yoon et al., again consisting of a rotat- 
ing core, hot envelope (which they term "corona"), and 
a Keplerian disk. Perhaps owing to different initial con- 
ditions, the final core temperatures differed by orders of 
mag nitude between the two groups. Loren-Aguilar et al. 



of this concern, the review of Hillebrandt and Niemeyer 
Q dismissed the double degenerate scenario. Such a col- 



li^ . Because [124j concluded that the disks surrounding the remnant 



are very turbulent, and argued that this would lead to 
high accretion rates and, subsequently, off-center igni- 
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tion and accretion induced collapse. Shen et al. [121 
also find that C does not explosively ignite during the 
merger, but for different physical reasons owing to the 
ability of magnetic stresses to redistribute angular mo- 
mentum. Over a longer timescale after the merger, they 
find that conditions at the interface of the tidally dis- 
rupted secondary and primary white dwarf are such that 
convective C burning should ensue leading to slow shell 
burning that eventually results in a high-mass O-Ne-Mg 
white dwarf or accretion induced collapse. 

Other work serves to stress the importance of the fi- 
delity of nu meri cal methods and in clud ed physics. Fryer 
and Diehl [l26l ] and Diehl et al. jl27l | performed sim- 
ulations of merging white dwarfs with similar methods 
but fo und differ ent re sults. The collaboration of D'Souza 
et al. [l28j and |l29l | with different, grid-based methods 
found that, depending on the initial configuration, sta- 
ble mass transfer c ould be achieved, a result s een o nly 
by Fryer and Diehl [Hi]. Recently, Dan et al. [HI in- 
vestigated the issue of stable mass transfer and found 
it possible wit h a c areful treatment of initial conditions 
and Dan et al. |l3f| performed a wide-ranging survey, but 
found the chance of a successful detonation following a 
merger of C-0 white dwarfs small. A nice opus on m any 
of the concerns was produced by Pakmor et al. [l32j . 

Recent studies indicate su ccess with white dwarf 
merger models. Pakmor et al. [l33j found that the case 
of nearly equal masses of w 0. 9M C .) may explain sub- 
luminous events. Pakmor et al. [116] found that violent 
mergers with primary masses of ~ 0.9Mq and mass ratios 
< 0.8 will produce explosions and are promising candi- 
dates for s ub-luminous events. Very recently, Pakmor 
et al. [l34| found that under the right conditions, the 
violent merger of two C-0 white dwarfs may produce a 
normal Type la supernova. 

The merger of two white dwarfs is also thought to 
poss ibly produce a rapidly rotating massive white dwarf 
[135| . Following considerable exploration of detona- 
tion scenarios including one-dimensional rotating and 
off-center models (sec Rcf. [72[ and references therein), 
the first multidimensional simulations of detonations in 
rapidly rota ting white dwarfs were performed by Stein- 
metz et al. [135J. The result of these two-dimensional 
simulations indicated that as with earlier studies of det- 
onations [txj, the bulk of the star burned all the way 
to Fe-group elements and did not produce the expected 
intermediate-mass elements. 

Recently Pfannes et al. jl36l | performed full three- 
dimensional simulations of turbulent deflagrations in 
rapidly rotating massive white dwarfs. They found 
that the amount of burning in their rotating models 
was similar to that of deflagrations occurring in non- 
rotating massive white dwarfs, producing "comparably 
weak and anisotropic explosions that leave behind un- 
burnt material at the center and in the equatorial plane." 
The conclusion drawn from their deflagration models is 
that rotation will not explain the obs erved variation in 
peak luminosities. Pfannes et al. 1371 ] explored detona- 



tions occurring in rapidly rotating massive white dwarfs 
with three-dimensional simulations. In this case, mod- 
els indicated that prompt detonations could produce 
bright supernova-like events, including producing suffi- 
cient amounts of intermediate-mass elements. They con- 
clude that detonations in rapidly rotating massive white 
dwarfs can explain some bright events. 

Additional support for the idea that rapidly rotating 
massive white dwarfs might explain some b right events 
comes from the recent study of Hachisu et al. |f38| . They 
presented a model that includes optically thick winds 
from an accreting white dwarf, mass stripping from the 
binary companion by the white dwarf wind, and support 
from differential rotation. They report that masses can 
reach up to 2.7M© and explored the results of explosions 
occurring in three mass ranges that depend on the char- 
acter of the rotation and the onset of secular instability 
and conclude that the model can explain both bright and 
normal events, and that differences in angular momen- 
tum loss rates (and thus spindown rates) might explain 
"prompt" and "tardy" components of supernova obser- 
vations. 



3. Double Detonation Models 

Another progenitor channel under study is known as 
the "double detonation" model. This idea is similar to 
the single degenerate scenario described in Section HP II 
but in this case the white dwarf does not need to ap- 
proach the Chandrasekhar limiting mass. Pioneering 
work introducing double d eton ation m odels was per- 
formed by Woosley et al. [l39j |. Taam [HJ [l4l| and 
Nomoto |l42l . Il43| . Studies indicated that the under 
the right conditions, an accreted He layer could flash 
with such strength as to produce an inwardly propa- 
gating shock that would ignite a detonation in the C- 
O core. Studies also indicated that the double deto- 
nation scenario could work for a wide range of white 
dwarf masses, not just the near- Chandrasekhar case 



144j | . and, accordingl y, low er mass models were dubbed 



"sub- Chandrasekhar" [145| 

Multidimensional simulations confirmed many of these 
findings, but as with other models, uncertainties in parts 
of the problem like the initial conditions were found to 
play a significant role on the explosion outcome. Ac- 
cordingly, a common conclusion was that do uble deto- 
nation models may explain some events [HI, Il46l - ll51 |. 
Most models included a massive He layer ~ 0.2Af© 
so that intermediate and heavy elements synthesized in 
the He detonation of such models would appear in the 
outer parts of the eject a, a result at odds with observa- 
tions P. ll48l . fl52lll53t . 

Recently, however, Bildstcn et al. 154| found that 
fairly thin He layers, ~ O.O6M , could become dynam- 
ically unstable and flash on sub - Chandrasekhar mass 
white dwarfs, and Fink et al. jl55| found that a He det- 
onation at the base of the accreted layer "robustly trig- 
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gcrs" seco ndary detonations in the C-0 core. Kromer 
et al . |l56l | further investigated the models of Fink et al. 
jl55l | and calculated light curves. They found the right 
range of rise and decline times and brightness, but found 
the color too red due t o com position of the He shell ash. 

Woosley and Kasen [157| | performed an extensive one- 
dimensional parameter study, calculating nucleosynthe- 
sis and light curves for a variety of white dwarf masses 
and temperatures and He accretion rates. They found 
that only hot, massive white dwarfs with thin He layers 
produce la-like events and discussed a number of poten- 
tial objections to these models as a general solu tion to 
the Type la supernova problem. Sim et al. |l53| inves- 
tigated central detonations in bare (no He layer) sub- 
Chandrasekhar mass C-0 white dwarfs. They produced 
light curves and detailed nucleosynthesis results and con- 
cluded that such models are "in qualitatively good agree- 
ment with observed properties of Type la supernovae," 
but note that the question of realistic progenitor syst ems 
in keeping with these models remains. Sim et al. [1581 ] in- 
vestigated double detonation scenarios with both edge- lit 
and central C detonations. They found that for relatively 
higher mass accreted He shells, conditions are likely for 
a core detonation to occur following a He detonation for 
all expected C-0 core masses. The light curves of these 
core detonation models were readily distinguishable from 
edge-lit or He detonation with no core detonation mod- 
els. All these recent studies serve to stress the importance 
of details of the accreted layer to the distribution of nu- 
clei in the ejecta and the need for high-fidelity spectral 
observations for validating this class of models. 



II. SINGLE DEGENERATE DELAYED 
DETONATION MODELS 



The single degenerate scenario posits comprcssional 
heating of a white dwarf due to accretion. The core of 
the white dwarf begins ru naway he ating due to C fusion 
outpacing neut rino losses |l59lll60j and a period of con- 
vection ensues [l6ll and references therein]. When heat 
generation outpaces transport of heat by convection, a 
flame is born that eventually will incinerate the white 
dwarf. The central density at flame igni tion is slightly 
decreased (see e.g. Piro and Bildsten Il62l ) and is gener- 
ally near ~2x 10 9 g cm -3 for successful models of Type 
la supernovae (see e.g. Refs. [52l |94|). 

A variation of the delayed detonation mechanism 
(within the single degenerate paradigm) is that of a de- 
flagration to detonation transition (DDT). After igni- 
tion, the flame propagates as a subsonic deflagration for 
a while and then transitions to a supersonic de tonation 
that rapidly consumes the star (921494 , Il0ll4l06| . We de- 
scribe the details of our implementation of this explosion 
mechanism below. 

At present, the physical mechanism by which a DDT 
in degenerate supernova material occurs is an area of cur- 
rent research (see Refs. 16314168 [ and references therein) . 



Simulations of supernovae involving DDT assu me it oc- 
curs via the Z eldovich- gradient mechanism (Ref . Il03l but 
see also Refs. Il0(j Il64 ). in which a gradient in reactivity 
leads to a series of explosions that are in phase with the 
velocity of a steadily propagating detonation wave. Many 
authors suggest that when the flame reaches a state of 
distributed burning, which is when turbulence on scales 
at or below the laminar flame width is fast enough to 
dominate transport processes (see e.g., |l69j |). fuel and 
ash are mixed and the temperature of the fuel is raised 
and "prepared" in such a way to produce the required 
reactivity gradient. A requirement for distributed burn- 
ing is that the ratio of turbulent intensity to the lami- 
nar flame speed must exceed so me unknown thres hold, 
which is still actively researched fiol Hoi [r70l - fT75j . En- 
trance into the distributed burning regime does not guar- 
antee such a react ivity gradient to form. Woosley [lT6| 
and Woosley et al. [165| studied incorporating more strin- 
gent requirements for these conditions to be met. 

For burning under conditions found in the white dwarf, 
the ratio of turbulent intensity to laminar flame speed 
changes most rapidly due to the change in laminar flame 
properties, which are strongly dependent on fuel den- 
sity. Therefore, a DDT is typically assumed to occur 
at a range of densities tha t generall y lie in the range of 
10 6 - 7 to 10 77 g cm" 3 [l03 . fl7^4l8H . although some re- 
search implements DDT criteria that include the strength 
of background turbulence (28l . and references therein] . In 
simulations, the choice of pddt increases or decreases the 
duration of the deflagration phase, which increases or de- 
creases the amount of expansion prior to the detonation 
and hence the co mpo sition of the yield of Iron Group 
Elements (IGEs) fl8ll ]. 



III. CONTEMPORARY RESEARCH IN THE 
DDT PARADIGM 

As outlined above, DDT models agree well with ob- 
servations and accordingly have been accepted as the fa- 
vored model of Type la supernovae. As the models work 
well, contemporary research goes well beyond questions 
of whether or not explosions happen to explore secondary 
parameters beyond the pri mary lig ht curve parameter, 
the brightness or Ami^{B) [271. [281] . In this section, we 
address the role of composition and central density on the 
explosion of models of Type la supernovae in the single 
degenerate DDT paradigm. 



A. The Role of Progenitor Composition 

The composition of a white dwarf can influence the 
explosion in a host of ways, including the density at 
which the deflagration ignites, the flame speed, the den- 
sity at which the DDT occurs, the energy released dur- 
ing the phases of burning, and the neutron excess of the 
burned material, which strongly influences properties of 
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the light curve of an event, including the peak bright- 
ness, expansion velocities, and the width-luminosity rela- 
tion. The most important constituents of the progenitor 
white dwarf are 12 C, 16 0, and 22 Ne. The composition is 
principally set by the post-main sequence evolution [29| . 
with the inner core formed during the convective core 
He-burning phase, and the outer layers f orme d in shell 
burning on the asymptotic giant branch [l82l and ref- 
erences therein]. These two burning stages operate in 
very different burning regimes with core He burning oc- 
curring at lower temperatures and higher densities than 
shell burning. At these lower temperatures in the core, 
the 12 C (a, 7) 16 reaction is favored over 3a, and as a 
result 12 C is depleted in the core, while it is not in the 
outer layers [Hj]. 

Additionally, the zero-age main sequence (ZAMS) 
mass and mctallicity influence the composition. As the 
ZAMS mass is increased from the low end of the range 
that produce C-0 white dwarfs, the central tempera- 
ture is increased during core He burning to favor 3a over 
12 C (a, 7) 16 0, thereby increasing the C/O ratio. As the 
ZAMS mass continues to increase, however, convection in 
the core becomes a dominant process in transporting heat 
away from the core and the competition between 3a and 
12 C (a, 7) 16 shifts back to favor 12 C (a, 7) 16 0, result- 
ing in a reduced C/O ratio for the highest ZAMS masses 
capable of producing CO white dwarfs. The C/O ratio is 
sensitive to the treatment of the conv ection process and 
is still actively researched [29|, [H, Il82| | . Increased metal- 
licity increases the opacity and results in a smaller core 
mass and lower central temperatu re, t hus mimicking the 
effects of a decreased ZAMS mass jl83l | . For a given tem- 
perature, increased metallicity also reduces the 3a rate 
to favor 12 C destruction via 12 C (a, 7) 16 0, which results 
in a lower C/O ratio; albeit, the effect is negligible for 
Z < 0.02 Hi. 

Neutron rich isotopes are synthesized in two epochs. 
First, during He burning, the aboriginal C, N, and O 
is transformed into 22 Ne HU, leading to a direct de- 
pendence on metallicity of the parent stellar population. 
In addition, neutron-rich material is formed in the pre- 
cxplosion c onvective C burning core at a comparable 
abundance fl62l.[l83. 



The density at which the ignition of the deflagration 
occurs is sensitive to several factors, including 12 C abun- 
dance in the core and the thermal history of the core, 
which follows from the accretion rate and possibly prop- 
erties of the He flashes (52[. Either of these may have 
metallicity dependencies involved with both the evolu- 
tion of the parent s tella r population and the still incom- 
pletely understood 185] process of progenitor system for- 
mation. There are also uncertainties in the screening 
enh ancement of nuclear reactions at these high densi- 
ties [MEH, 

and in the reaction rates themselves, par- 
ticularly the 12 C + 12 C reaction rate [77], and references 
therein] . 

The energy per unit mass released during thermonu- 
clear burning influences both the dynamics of the white 



dwarf expansion during the deflagration phase and the 
dynamics of local buoyancy effects that accelerate the 
flame during the deflagration phase. The composition 
influences the energy release of thermonuclear burning 
in two ways. First, the gross nuclear energy available 
per unit mass is mostly sensitive to the abundance of 
12 C. Second, by changing the balance of protons and 
neutrons, the 22 Ne abundance influences the "ash" to 
which the material burns, with neutron rich material fa- 
voring more tight ly bound nuclei and therefore releasing 
more energy jl88l ] . The abundances over a broad region 
of the white dwarf are involved in this case, and, accord- 
ingly, involve the products of both the He core and shell 
burning. 

The deflagration propagates via thermal conduction 
and the laminar flame speed is sensitive to both the 12 C 
and 22 Ne abundances that influence the energy release 
and t he s peed of the early stages of the nuclear fusion 
[l89l Il90| . Turbulence from the simmering phase and 
resulting from buoyancy instability as the flame propa- 
gates has the ability to add a nearly arbitrary amount of 
local area to the flame surface. The effect is to make the 
burning rat e lar gely independent of the laminar propaga- 
tion speed [l9l| . The early deflagration phase, however, 
when the laminar flame speed interacts with the lower 
strength turbulence from core convection sets the mor- 
phology of the burned regi on a t the point when strong 
buoyancy effects take over [l92| . 

Although the physics of DDT are incompletely under- 
stood, the process will depend on composition because 
both the flame speed and width depend on the abun- 
dances of 12 C and 22 Ne. Depending on the details of 
the DDT, the transition could occur either in material 
with the core composition or with the outer layer com- 
position due to the nature of subsonic burning and the 
response of the star. Figure [T] shows four snapshots of 
a two-dimensional simulation in which DDT is triggered 
at a specified density of 10 71 gcm -3 (green contour) in 
a progenitor white dwarf with a different composition in 
the core than the outer layer (marked by the blue con- 
tour). As the deflagration ensues, the star expands and 
the density contour shrinks. For this transition density, 
DDT occurs in material from the core, but for lower tran- 
sition densities, it may be possible for DDT to occur in 
material from the outer layer. The core 12 C abundance 
will be lower than the outer layer due to the differences in 
temperature experienced in shell burning as opposed to 
core He burning. The neutron excess in the core is also 
expected to be greater than the outer layer due to the 
slow C burn i ng phase as w ell as potential sedimentation 
effects fl62l HM [193I4T95] . Because the DDT density 
determines the duration of the deflagration phase and 
thus degree of expansion of the star at the detonation, 
the influence of metallicity on the DDT densi ty p l ays a 
significant role in the outcome of the explosion |l8ll . Il84| . 

In addition to the indirect effect on energy release de- 
scribed above, the 22 Ne content, which sets the neu- 
tron excess, has a direct influence over the final nucle- 



9 



osynthetic products, particularly the amount of 56 Ni. 
Timmes et al. [BTJ showed that the decrease in the 56 Ni 
produced in the explosion, absent other factors, should 
be linearly proportional to the 22 Ne content and there- 
fore the original metallicity of the stellar po pulat ion, a 
result seen in some multidimensional results 196]. The 
distribution of 22 Ne within the white dwarf is important, 
with 22 Ne in the interior of the star influencing the gross 
yield and 22 Nc in surface layers influencing the ejecta 
opacity. Neutron excess also influences the structure of 
a white dwarf because it is supported by the pressure of 
degenerate electrons. The neutron excess, determined by 
the 22 Ne abundance, sets the weight of baryons each elec- 
tron must support. Accordingly, white dwarfs of a given 
mass with a higher neutron excess will be more compact 
because there are fewer total electrons. Conversely, at 
the same central density, a star with more 22 Ne and thus 
more neutrons, will be slightly less massive. These con- 
current events on th e den sity and structure of the white 
dwarf may be small |l97l | , but should be included in re- 
alistic models due to the marginally bound nature of a 
near-Chandrasekhar mass white dwarf. 

In a comprehensive st udy using one dimensional DDT 
models, Hocflich et al. |l05j addressed the influence of 
initial composition on the nucleosynthesis, light curves, 
and spectra along with cosmological implications. They 
found a decreasing C/O ratio released less energy during 
burning, which increased the duration of the deflagra- 
tion phase, which in turn decreased the production of 
56 Ni and influenced both the bolometric and monochro- 
matic light curves. They also varied the DDT transition 
density in one model, and found that reduction of the 
transition density produced similar results on the yield 
to decreasing the C/O ratio. They reported, however, 
that there arc differences in the distribution of Si in the 
expanding remnant that lea d to s light differences in the 
light curves. Hocflich et al. f 105| also reported that in- 
creasing the metallicity of the progenitor decreased the 
yield of 56 Ni because the decreased number of protons 
led to the production of relatively more stable Fe-group 
nuclei. 

W hile one-dimensional simulations by Umcda et al. 
jl98l | also found that a reduced C/O ratio decreased the 
production of 56 Ni, later three-dimensional defiagration- 
only simulations by Ropke and Hillebrandt jl99l | found 
that the C/O ratio does not significantly influence the 
final yield of 56 Ni. They attribute this to a complex 
interplay between the nucleosynthesis and nonlinear tur- 
bulent flame evolution. This result was later confirmed 
with higher resolution simulations and post-processing 
tracer particles with a more detailed nuclear reaction net- 
work [261 ]; however, due to the highly non-linear nature of 
the turbulent flame propagation, more simulations with 
different ignition conditions are necessary to rule out any 
systematic trend. The effect of the C/O ratio on DDT 
has y et to be explored in detail; although, Umcda et al. 
|198| speculated that a larger C/O ratio should increase 
the transition density, and hence, increase the yield of 



56 Ni. However, if DDT occurs due to the Zeldovich- 
gradient mechanism in the distributed burning regime, 
then the transition density should increase with decreas- 
ing C/O ratio becaus e of t he s t rong er effect on the lam- 
inar flame properties jl8ll 1 190L I20C 



The results of Hocflich et al. [1051 ] contain an impor- 
tant point concerning DDT models we wish to stress. 
The amount of expansion of the white dwarf during the 
deflagration phase significantly influences the explosion 
outcome, particularly the mass of Fe-group elements. We 
will quantify this result below in the discussion of our re- 
search into the the role of the DDT transition density. 

Hoflich et al. (27j addressed secondary parameters of 
Type la supernova light curves, meaning parameters 
other than Am^B), that influence the brightness of 
an event. They argued that high quality lig ht curves 
obtained by the Carnegie Supernova Project |20ll |202| 
show clear evidence of the existence of secondary pa- 
rameters responsible for variations in the light curves 
and explored the relationship between possible param- 
eters, central density, the C/O ratio (directly follow- 
ing from the main sequence mass), and metallicity, and 
the light curves from explosion models. The study was 
based on DDT models with a fixed transition density 
of 2.3 x 10 7 gcm~ 3 . They reported that central den- 
sity and the C/O ratio may be treated as independent 
secondary parameters characterizing the light curve, and 
that larger main sequence masses (lower C/O ratios) pro- 
duces a slower rise of the light curve. 



B. The Role of Central Density 

As described above, the composition of the progeni- 
tor white dwarf can play an important role in the out- 
come of an explosion. The composition is determined by 
many factors, such as the mass of the main sequence star 
that became the white dwarf, but it is also influenced by 
properties of the host galaxy such as metallicity. These 
properties lead to correlations between galaxy type and 
brightness of events. A similar property of the progeni- 
tor white dwarf that influences the brightness and follows 
from the age of the progenitor system, and hence the age 
of the galaxy, is the central density. 

As the accreting white dwarf approaches the Chan- 
drasekhar mass, it is relatively insensitive to variations in 
the central density because of the degenerate equation of 
state. The configuration of the convecting interior is un- 
certain, however, because of an incomplete understand- 
ing of the physics that influence ignition. T here hav e 

been many studies of the convective phase |l6llll62l . l203l - 
|208| , but examples of remaining uncertainties include the 
12 C- 12 C rate and its influence on the ignition density 



[77l [209l - [2TT t and the loss of energy du e to neutrinos, i.e. 
the convective Urea process 212U 2 



217 



While these many questions remain, it is accepted that 
the initial mass of the white dwarf (i.e., the mass at which 
it formed) and the accretion history, particularly the pc- 
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riod of cooling prior to the onset of accretion, largely de- 
termine the state of the interior of the white dwarf [28| . 
The study of Lesaffre et al. [3(| found that white dwarfs 
experiencing a long cooling period prior to the onset of 
accretion cool to a lower central temperature, so that 
once accretion begins to reheat the core a greater increase 
in density is necessary to reach the conditions for a ther- 
monuclear runaway (see especially Figure 3 of Ref . l30h . 
Therefore, while the central density of the white dwarf 
may have been the same prior to the cooling period, at 
ignition of the thermonuclear runaway the central density 
is higher for white dwarfs with a longer cooling time. 

The effect of density on electron capture rates has been 
known for some time [90T and references therein]. Elec- 
tron capture rates increase with density, so that at the 
higher densities found in near-Chandrasekhar mass white 
dwarfs, the rates are rapid enough that significant neu- 
tronization occurs during the explosion. The issue of the 
effect of central densi ty has be en explored by many re- 
search groups [H-il, EJ, Elj| and below wc highlight 
some of this work. 



Hoflich et al. |27| in their study of secondary param- 
eters of Type la supernova light curves considered cen- 
tral density. They found that the central density influ- 
ences the later evolution of the light curve (~ 30 days 
after maximum and later) with increases in central den- 
sity causing the later portion of the light curve to shift 
down and vice versa. They attribute this difference to 
differences in the amount and distribution of 56 Ni in the 
central envelope. 

Recently Seitenzahl et al. [28| explored the effect of 
changing central density on the explosion outcome with 
three-dimensional models and a treatment of the energet- 
ics of the thermonuclear burning similar to the method 
wc employ (described below) that includes weak interac- 
tions, specifically electron capture, allowing their models 
to neutronize. Seitenzahl et al. [28| found that the in- 
creased neutronization rates decreased the relative abun- 
dance of 56 Ni in the Fe-group material synthesized in 
the explosion. Considering an identical spatial distri- 
bution of ignition kernels, however, they found the in- 
creased gravitational acceleration following from a higher 
central density increases the production of turbulence, 
which produces a higher burning rate and triggers the 
DDT sooner, with less expansion of the star and there- 
fore a higher yield of Fe-group elements. These two ef- 
fects combine to give the relatively constant yield of 56 Ni, 
and they conclude that the net effect is that the mass of 
56 Ni synthesized in an explosion appears insensitive to 
the central density of the progenitor white dwarf. More 
recent three-dimensional results including detailed nucle- 
osynthetic post-processing of 1 x 10 6 tracer particles per 
model suggest that 56 Ni mass increases with central den- 
sity. In those models, the fraction of 56 Ni of IGEs de- 
creases due to neutronization, but the increase in the 
total amou nt in IGE outweighs the decrease due to neu- 
tronization [220j ] . 

Seitenzahl et al. [28[ also report that the 56 Ni yield 



is sensitive to the initial conditions, which largely de- 
termine the behavior of the deflagration phase. They 
found that a strong deflagration phase following from 
a many-kernel (larger) ignition region produces a lower 
yield of 56 Ni primarily because of the increased expansion 
of the star during the deflagration phase but also due to 
increased neutronization during the deflagration phase. 
They note that this confirms previous studies indicating 
that the strength of the defla gratio n phase is a primary 
parameter of 56 Ni production |22l| . 



IV. RESEARCH RESULTS 

In this section we provide a brief summary of our own 
theoretical research into systematic effects on the bright- 
ness of Type la supernova explosions. We describe our 
methodology, including our models and statistical ap- 
proach to studying systematic effects, and present high- 
lights from our investigation into the role of the metal- 
licity and central density of the progenitor white dwarf. 
We gratefully acknowledge contributions to the results 
presented here from our many collaborators. 



A. Methodology 

Our methodology for the theoretical study of Type 
la supernovae consists of four principal parts. First 
is the ability to construct parameterized, hydrostatic 
initial white dwarf progenitors that can freely change 
thermal and compositional structure to mat ch fe atures 
from the literature about progenitor models 181]. Sec- 
ond is a model flame and energetics scheme with 
which to track both (subsonic) deflagrations and (su- 
personic) detonations as well as the evolution of dy- 
namic ash in nuclear statistical equilibrium (NSE). This 
flame/energetics sc heme is im plemented in the Flash hy- 
drodynamics code [222h224| . Third is utilization of a 
scheme to post-process the density and temperature his- 
tories of Lagrangian tracer particles with a detailed nu- 
clear networ k in orde r to calculate detailed nucleosyn- 
thetic yields [225l |226| . Fourth, we developed a statisti- 
cal framework with which to perform ensembles of sim- 
ulations for well-controlled studies of systematic effects. 
Below we highlight the flame model and the statistical 
framework. Complete details of the meth odology can b e 



found in previously published results [H [HI [225M227 



The simulations presented here are all two-dimensional. 



1. Flame Model 

The significant disparity between the length scale of a 
white dwarf (~ 10 9 cm) and the width of laminar nu- 
clear flame (< 1 cm) necessitates the use of a model 
flame in simulations of Type la supernovae. Even sim- 
ulations with adaptive mesh refinement cannot resolve 
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the actual diffusive flame front in a simulation of the 
event. The model we use propagates an artificially broad- 
ened flame fro nt w i th a n advection-diffusion-rcaction 
(ADR) scheme [228l |229| | that has been demonstrated 
to be aco ustically quiet and produce a unique flame 
speed |227| . This scheme evolves a reaction progress vari- 
able 4>, where = indicates unburncd fuel and <f> = 1 in- 
dicates burned ash, with the advection-diffusion-rcaction 
equation 

d t (j) + u- V0 = kV 2 + -R{4>) ■ (1) 

T 

Here k is a constant and R((f>) a non-dimensional func- 
tion, and both are tuned to propagate t he r e actio n front 
at the physical speed of the real flame [l90l |23C| and to 
be just wide enough to be resolved in our simulation. We 
use a modified versio n of t he KPP reaction rate discussed 
by Vladimirova et al. |229| , in which R oc (</>— e)(l — 0+e), 
where e ~ 10~ 3 , w hich is acoustically quiet and gives a 
unique flame speed [2271 ]. 

In simulating the deflagration phase, the flame front 
separates expanded burned material (the hot ash) from 
denser unburncd stellar material (cold fuel). The ex- 
pansion and buoyancy of the burned material forces the 
interface upward into the denser fuel, and the configura- 
tion is susceptible to the Rayleigh- Taylor fluid instabil- 
ity [23ll . |232| . It is necessary to enhance the flame speed 
in order to prevent turbulence generated by the Raleigh- 
Taylor instability from destroying the artificially broad- 
ened flame front. In the simulations discussed here, the 
enhancement is ac complished by the method suggested 
by Khokhlov [228| | in which we prevent the flame front 
speed, s, from falling below a threshold that is scaled with 
the strength of the Raleigh- Taylor instability on the scale 
of the flame front (s oc y/gl, where g is the local gravity 
and and £ is the width of the artificial flame, which is a 
few times the grid resolution). This scaling of the flame 
speed prevents the Rayleigh- Taylor instability from ef- 
fectively pulling the flame apart and also mimics what 
is a real enhancement of the burning area that is occur- 
ring due to structure in the flame surface on unresolved 
scales. It is expected that, for much of the white dwarf 
deflagration, the flame is "self-regulated", in which the 
small scale structure of the flame surface is always suf- 
ficient to keep up with the large-scale buoyancy-driven 
motions of the burned material. Thus the actual burn- 
ing rate is determi ned by this action, which is resolved 
in our simulations |233l |. 

This technique has demonstrated a suitab le level o f 
convergence for studies of Type la supernovae |227l |234| . 
The technique explicitly drops the flame speed to zero 
below a density of 10 7 g cm~ 3 , approximately where the 
real flame will be extinguished. This prescription cap- 
tures some effects of the Rayleigh- Taylor instability and 
maintains the integrity of the thickened flame, but it does 
not completely describe the flame-turbulence interaction. 
Also, we neglect any enhancement from background tur- 
bulence from convection prior to the birth of the flame. 



Future work will i nclude physically-motivated models for 
these effects j235l |. 

In addition to the model flame, the energetics scheme 
also includes the nuclear energy release occurring at the 
flame front, in subsequent burning stages, and in the dy- 
namic ash in nuclear statistical equilibrium. We per- 
formed a detailed study of the nuclear processes occur- 
ring in a flame in the interior of a white dwarf and de- 
veloped an efficient and accurate method for incorporat- 
ing the results into numerical simulations jU Il88l |227| . 
Tracking even tens of nuclear species is computationally 
prohibitive, and many more than this are required to ac- 
curately calculate the physics such as electron capture 
rates that are essential to studying rates of neutroniza- 
tion. We instead reproduce the energy release of the 
nuclear reactions with a highly abstracted model based 
on tabulation of properties of the burned material calcu- 
lated in our study of the relevant nuclear processes. This 
method accurately captures the thermal history of the 
material as it burns and evolves, which enables embed- 
ded particles to obtain accurate Lagrangian density and 
temperature histories. Detailed abundances can then be 
recovered by post-processing these time histories with a 
nuclear network including hundreds of nuclides. 

The nuclear processing can be well-approximated as a 
three stage process: initially C is consumed, followed by 
O, which creates a mixture of Si group and light elements 
that is in quasi-statistical equilibrium, also known as 
nuclear statistical quasi-c quilibrium (NSQE) [236M240| 1 : 
finally the Si-group nuclei are converted to Fc group, 
reaching full nuclear statistical equilibrium, NSE. In both 
of these equilibrium states, the capture and creation of 
light elements (via photodisintegration) is balanced, so 
that energy release can continue by changing the rela- 
tive abundance of light (low nuclear binding energy) and 
heavy (high nuclear binding energy) nuclides, an action 
that releases energy as buoyant burned material rises and 
expands. We track each of these stages with separate 
progress variables and separate relaxati on t imes derived 
from full nuclear network calculations [l88j . We define 
three progress variables representing consumption of C, 
4>c consumption of O to material in NSQE, 0nsqe, and 
conversion of Si to Fe-group nuclides to form true NSE 
material, 0nse- The physical state of the fluid is tracked 
with the electron number per baryon, Y e , the number of 
nuclei per baryon, Y\, and the average binding energy per 
baryon, q, the minimum properties necessary to hydro- 
dynamically evolve the fluid. C consumption is coupled 
directly to the flame progress variable, 0c = 0j from eq. 
(H|) above, and the later flame stages follow using simple 
relationships from more detailed calculations. 

Burning and evolution of post-flame material change 
the nuclear binding energy, and we use the binding energy 
of magnesium to approximate the intermediate burning 
products of C. The method is finite differenced in such a 
way to ensure explicit conservation of energy. Weak pro- 
cesses that neutronize the material (e.g. electron capture) 
are included in the calculation of the energy input rate, 
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as are neutrino losses, which are calculated by convolv- 
ing the NSE distribution with the weak interaction cross 
sections. Tracking the conversion of Si-group nuclides 
to the Fe-group is important for studying the effects of 
electron capture because the thresholds are lower for the 
Fe-group nuclides. Both the NSE state and the electron 
capture rates were calculated with a set of 443 nuclides 
including all which have weak interactio n cro ss sections 
given by Langankc and Martfnez-Pinedo |24l| . 

Electron capture feeds back on the hydrodynamics in 
three ways: the NSE can shift to more tightly bound 
elements as the electron fraction, Y e , changes, releasing 
some energy and changing the local temperature; also 
the reduction in Y e lowers the Fermi energy, reducing 
the primary pressure support of this highly degenerate 
material and having an impact on the buoyancy of the 
neutronized material; finally neutrinos are emitted (since 
the star is transparent to them) so that some energy is 
lost from the system. Also, as we will see below, increased 
rates of neutronization produce yields of more neutron 
rich material and therefore less 56 Ni, thereby influencing 
the brightness of an event. Complete details of the NSE 
calculations may be found in Ref. Seitenzahl et al. [242j 
and the details of the implementation in our simul ation s 
may be found i n Tow nsley et al. [3l| , Calder et al. |l88| , 
Townsley et al. (22^ . 

In addition to all of these effects during the deflagra- 
tion phase of Type la supernovae, the progress- variable- 
based method has bee n extende d to model the gross 
features of detonations j226l |243| . Instead of coupling 
the first burning stage, <pc , representing C consumption, 
to the ADR flame front, we instead can use the actual 
temperature-dependent rate for C burning, or possibly a 
more appropriate effective rate. Doing so allows shock 
propagation to trigger burning and therefore create a 
propagating deton ation . This method has been used suc- 



cessfully by Ref. |244j in modern studies of the DDT, 
and our multistage burnin g m odel s hare s many features 
with theirs (see also Refs. 92 1 and [245j ) . We treat the 



later burning stages very similarly, though we have taken 
slightly more care to track the intermediate stages and 
have nearly eliminated acoustic noise when coupling en- 
ergy release to the flame. 



2. Statistical Framework 

We also developed a theoretical framework for the 
study of systematic effects in Type la supernovae [3l| . 
This framework allows the evaluation of the average de- 
pendence of the properties of supernovae on underlying 
parameters, such as composition, by constructing a the- 
oretical sample based on a probabilistic initial ignition 
condition. Such sample-averaged dependencies are im- 
portant for understanding how Type la models may ex- 
plain features of the observed sample, particularly sam- 
ples generated by large dark energy surveys utilizing 
Type la supernovae as distance indicators. 



The theoretical sample is constructed to represent sta- 
tistical properties of the observed sample of Type la su- 
pernovae such as the mean inferred 56 Ni yield and vari- 
ance. Within the DDT paradigm, the variance in 56 Ni 
yields can be explained by the development of fluid in- 
stabilities during the deflagration phase of the explosion. 
By choice of the initial configuration of the flame, we may 
influence the growth of these fluid instabilities resulting 
in varying amounts of 56 Ni synthesized during the explo- 
sion. In Townsley et al. (3l| , we found that perturbing a 
spherical flame surface with radius (tq = 150 km) with 
spherical harmonic modes (Y™) between 12 < I < 16 
with random amplitudes (A) normally distributed be- 
tween — 15 km and, for three-dimensional simulations, 
random phases (<5) uniformly distributed between -ir — n 
best characterized the mean inferred 56 Ni yield and sam- 
ple variance from observations: 



id) =r + £ AiYiiO) 



(2) 



(max I 



ATe<' 



(=( m in m——l 



With a suitable random-number generator, a sample 
population of progenitor white dwarfs is constructed by 
defining the initial flame surface for a particular progen- 
itor. 



B. Investigation into Metallicity 

We investigated the direct effect of reprocessed stellar 
material (metals) in the host galaxy via the initial neu- 
tron excess of the progenitor white dwarf in Townsley 
et al. [31| . Because of weak interactions, metals pro- 
duced by nuclear burning are more neutron rich than H 
and He, and, accordingly, the neutron excess of these 
elements is thought to drive the explosion yield toward 
stable Fe-group elements. Thus, there is relatively less 
radioactive 56 Ni in the NSE mix, which results in a dim- 
mer event [5l|. We investigated this effect by introduc- 
ing 22 Ne into the progenitor white dwarf as a proxy for 
(neutron-rich) metals. The presence of 22 Ne influences 
the progenitor structure, the energy release of the burn, 
and the flame speed. The study was designed to mea- 
sure how the 22 Ne content influences the competition be- 
tween rising plumes and the expansion of the star, which 
determines the yield. We performed a suite of 20 DDT 
simulations varying only the initial 22 Ne in a progenitor 
model, and found a negligible effect on the pre-detonation 
expansion of the star and thus the yield of NSE mate- 
rial. The neutron excess sets the amount of material in 
NSE that favors stable Fe-group elements over radioac- 
tive 56 Ni. Our results were consistent with earlier work 
calculating the direct modification of 56 Ni mass from ini- 
tial neutron excess 1511. 



13 



In Jackson et al. 181| we expanded the earlier 
study [3l| to include the indirect effect of metallicity in 
the form of the 22 Ne mass fraction through its influence 
on the density at which the DDT takes place. We sim- 
ulated 30 realizations each at 5 transition densities be- 
tween 1 — 3 x 10 7 g cm~ 3 for a total of 150 simulations. 
We found a quadratic dependence of the NSE yield on 
the log of the transition density, which is determined by 
the competition between rising unstable plumes and stel- 
lar expansion. By then considering the effect of metal- 
licity on the transition density, we found the NSE yield 
decreases slightly with metallicity, but that the ratio of 
the 56 Ni yield to the overall NSE yield does not change 
as significantly. Observations testing the dependence of 
the yield on metallicity remain somewhat ambiguous, but 
the d epen dence we found is comparable to that inferred 
from [200]. We also found that the scatter in the results 
increases with decreasing transition density, and we at- 
tribute this increase in scatter to the nonlinear behavior 
of the unstable rising plumes. 

Figure [2] illustrates these results by plotting mass of 
56 Ni vs. metallicity for our m odels and the observational 
results of Konishi et al. 246] . The constant DDT density 
results (blue curve) represent the direct effect of intro- 
ducing 22 Ne into the progenitor white dwarf as a proxy 
for (neutron-rich) metals. Th e red curves present results 
adapted from Jackson et al. |l8l| that include the indi- 
rect effect of metallicity influencing DDT density. The 
results inferred from Konishi et al. [240] arc the black 
points. In the plot, metallicity is quantified by compar- 
ing the ratios of Fe to H in the progenitor white dwarfs 
relative to our Sun |247{ . We write this relationship as 



[Fe/H] = log 



10 



(AW^h) Sl 



C. Investigation into Central Density 

We investigated the effect of cent ral density on the ex- 
plosion yield in Krueger et al. 248j . While we found lit- 
tle change in the overall production of NSE material, we 
found a definite trend of decreasing 56 Ni production with 
increasing progenitor central density. We attribute this 
result to higher rates of weak interactions (electron cap- 
tures) that produce a higher proportion of neutronized 
material at higher density. More neutronization means 
less symmetric nuclei like 56 Ni, and, accordingly, a dim- 
mer event. This result may explain the observed decrease 
in Type la supernova brightness with increasing delay 
time. The central density of the progenitor is determined 
by its evolution, including the transfer of mass from the 
companion. If there is a long period of cooling before the 
onset of mass transfer, the central density of the progen- 
itor will be higher when the core reaches the C ignition 
temperature [30j, thereby producing less 56 Ni and thus a 
dimmer event. In addition, we found considerable vari- 
ation in the trends from some realizations, stressing the 



importance of statistical studies. 

In Krueger et al. [249j ] we expanded our earlier investi- 
gation into the role of central density on the production of 
56 Ni and hence the brightness of an event. We presented 
the details of our models, and found that the distribu- 
tion of 56 Ni in our models typically shows a somewhat 
clumpy concentration in the inner region of the remnant, 
with some dependence of the morphology on the cen- 
tral density. We found a trend of a shortened period of 
deflagration at higher central densities, but owing to in- 
creased rates of neutronization, the total neutronization 
is greater at higher densities, which explains our trend 
of decreasing 56 Ni production with increasing progenitor 
central density. We also considered main sequence evo- 
lution in comparing the masses of 56 Ni produced to ob- 
servations thereby obtaining an improved brightness-age 
relation from our results. 

Figure [3] presents results adapted from Krueger et al. 

j249| | relating the results to observations. The principal 
result of a simulation is the mass of 56 Ni produced, which 
is directly related to the brightness of an event. Using 
the method outlined in Howell et al. [HI, we can con- 
vert 56 Ni masse s to s tretch values, a measure related to 
peak brightness [250| , for a more direct comparison with 
observations. We combined the central density values 
with the results of Lesaffre et al. [3(j , which correlate the 
central density at the time of the ignition of the flame 
front to the cooling time of the progenitor WD (that is, 
the period of isolation prior to the onset of accretion), 
allowing us to express our results as ages. The results 
of Lesaffre et al. [30| suggest that a WD with a central 
density of 1 x 10 9 g cm" 3 will not ignite without further 
accretion, so for this comparison simulations with that 
value of central density were neglected. We also applied 
a shift to our data to account for the main sequence age, 
Tms> which we take to be constant across our results. We 
consider two estimated limiting values for tms (0.05 and 
1.0 Gyr, corresponding to main-sequence masses of ap- 
proximately 8 .0 a nd 1.5 solar masses, respectively; see 
Hansen et al. l25lh and included them with the original 
TMS = result in Figured] Adding in a tms consistent 
with our C-0 progenitors brings our results into better 
agreement with the two right-most points of Neill et al. 

Figure [3] shows that our theoretical results are not 
in complete agreement with the observed data. Obser- 
vationally, the age-brightness correlation may flatten at 
young ages, while our data do not, resulting in our data 
being too bright relative to young Type la supcrnovae. 
This study isolated the effects of central density and 
related that to age assuming a constant main-sequence 
mass, but there are other effects that may be correlated. 
Examples of such potentially correlated effects are: main- 
sequence mass and its correlation with central density, 
metallicity of the progenitor, core 12 C fraction prior to 
ignition of the deflagration, sedimentation, and others. 
Inclusion of such effects may modify the results presented 
here and are the subjects of future work. 
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V. SUMMARY AND CONCLUSIONS 

We have presented an overview of Type la supernova 
explosion models currently under study and a brief de- 
scription of contemporary research. Both observations 
and theoretical models are at a level of sophistication 
that allows comparison at a far more detailed level than 
just bulk properties such as energy release or bolometric 
light curve. Instead, studies explore secondary parame- 
ters of light curves and the role of properties such as the 
age and composition of the progenitor system on the ex- 
plosion yield and light curve. The structure and compo- 
sition of a white dwarf follow from properties of the host 
galaxy such as metallicity and age, and both theorists 
and observers are able to study trends of brightness with 
host galaxy type. These studies will better inform future 
observational campaigns probing the expansion history 
of the Universe and the properties of dark energy. 

Certainly, however, the problem is far from solved. The 
problems of discerning the progenitor systems (s) and ex- 
plaining the diversity of events remain. Researchers are 
actively exploring multiple progenitor systems, to explain 
"normal" events as well as outliers, both dim and bright. 
Modeling explosions in these progenitor systems is a chal- 
lenging problem, necessitating the development of com- 
plex sub-grid-scale models, and fundamental questions 
remain about models and the basic physics therein, e.g. 
the unsolved problem of turbulent combustion. Good 
progress is being made, however, and our assessment is 
that many of the outstanding questions will be answered 
in the near future. 

We presented results from our research into systematic 
effects on the production of 56 Ni and hence the brightness 
of an event. Our simulations assume the single degener- 
ate progenitor system and that the explosion proceeds via 
a DDT. Our approach to studying systematic effects is 
via suites of simulations allowing quantification of statis- 
tical properties. One important result from our research 
we stress here is the need for a statistical approach to dis- 
cerning trends as illustrated by the relatively large stan- 
dard deviation in Figure [2J 

Our models that include the direct effect of metallic- 
ity indicate a weak trend of decreasing brightness with 
increasing metallicity that is consistent with earlier re- 
sults [ll|. Models including the indirect effect of metal- 
licity on flame speeds and the deflagration to detonation 
transition density demonstrate a significantly stronger 
trend of decreasing brightness with increasing metallic- 
ity. Our results indicate a quadratic dependence of the 
yield of 56 Ni on the log of the transition density. Relat- 
ing metallicity to the transition den sity g ave relationship 



between brightness and metallicity j!81 



We also presented results from our research into the ef- 
fect of central density. Our suites of simulations showed a 



strong trend of decreasing brightness with increasing cen- 
tral density that we attribute to increase rates of neutron- 
ization. Relating our central density results to the age of 
the progenitor system allows comparison to observations, 
and in so doing we were able to provide a theoretical ex- 
planation of the observed trend that Type la supernovae 
from older host galaxies are systematically dimmer [249j . 
The agreement between our trend and the observational 
trends is not perfect, however, and we recognize that un- 
certainties abound and our two-dimensional models are 
incomplete. Plans for the future include targeted three- 
dimensional simulations including more complete physics 
to confirm our results. 

Our conclusion is that this is an exciting time to be 
studying the complex problem of thermonuclear super- 
novae. Substantial progress is being made by both ob- 
servers and theorists, and we hope that this review of one 
part of the problem, the role of chemical composition on 
explosion models, conveys some of this excitement. In 
preparing this review, we were struck by the remarkable 
progress that has been made in modeling this complex, 
multi-scale, multi-physics problem by many talented sci- 
entists. We also hope this review conveys our enthusiasm 
and respect for all of the careful study, both past and 
present, of the problem of Type la supernovae. 
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FIG. 1. A two-dimensional delayed-detonation simulation in 
which DDT is triggered at a specified density of 10 71 gcm -3 
shown by the green contour. The core of the progenitor white 
dwarf is more carbon-depleted than the outer layer with the 
core-outer layer boundary highlighted with the blue contour. 
Four different snapshots in time are provided that show the 
evolution of the flame and the expansion of the white dwarf. 
Note that the rightmost frame is scaled by a factor of two 
relative to the others. 
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Pddt = PDDT( 22Ne ) 
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FIG. 2. Plot of mass of Ni vs. metallicity comparing simula- 
tions that varied DDT densities with metallicity to the obser- 
vational results of Konishi et al. [246| (black points). The con- 
stant DDT result, blue curve, assumed that the DDT density 
did not vary with metallicity but included other metallicity 
effects on the structure of the progenitor star. The red curves 
present results from adapted from Jackson et al. |l8l| . with 
the the solid curve presenting the average value of simula- 
tions at a given metallicity and the widely-spaced dashed and 
dashed curves representing standard deviation and standard 
error, respectively. 
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FIG. 3. Plot of stretch vs. age comparing the scaled results 
of simulations varying the central density of the progenitor 
white dwarf to the observations of Neill et al. [i^] . In red are 
the points from this study with no shift (i.e., tms = Gyr), 
along with the standard error of the mean and a best-fit trend 
line. The green shaded region shows our best-fit line with 
tms = 0.05 — 1.0 Gyr. In blue are the binned and averaged 




